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ABSTRACT 

We investigate the effect of three important processes by which AGN-blown bubbles 
transport material: drift, wake transport and entrainment. The first of these, drift, 
occurs because a buoyant bubble pushes aside the adjacent material, giving rise to 
a net upward displacement of the fluid behind the bubble. For a spherical bubble, 
the mass of upwardly displaced material is roughly equal to half the mass displaced 
by the bubble, and should be ^ 10''"^ depending on the local ICM and bubble 
parameters. We show that in classical cool core clusters, the upward displacement by 
drift may be a key process in explaining the presence of filaments behind bubbles. A 
bubble also carries a parcel of material in a region at its rear, known as the wake. The 
mass of the wake is comparable to the drift mass and increases the average density of 
the bubble, trapping it closer to the cluster centre and reducing the amount of heating 
it can do during its ascent. Moreover, material dropping out of the wake will also 
contribute to the trailing filaments. Mass transport by the bubble wake can effectively 
prevent the build-up of cool material in the central galaxy, even if AGN heating does 
not balance ICM cooling. Finally, we consider entrainment, the process by which 
ambient material is incorporated into the bubble. Studies of observed bubbles show 
that they subtend an opening angle much larger than predicted by simple adiabatic 
expansion. We show that bubbles that entrain ambient material as they rise, will 
expand faster than the adiabatic prediction; however, the entrainment rate required 
to explain the observed opening angle is large enough that the density contrast between 
the bubble and its surroundings would disappear rapidly. We therefore conclude that 
entrainment is unlikely to be a dominant mass transport process. Additionally, this 
also suggests that the bubble surface is much more stable against instabilities that 
promote entrainment than expect for pure hydrodynamic bubbles. 
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1 INTRODUCTION 

The hot, gaseous atmospheres of galaxy clusters often show 
depressions in the X-ray surface brightne ss (see, for exam- 
ple McNamara 2002 ; Birzan et al. 2004; Dunn et al.1 l2005l : 



iRaffertv et al.l |2006| ) . These depressions are indicative of 
empty cavities, or bubbles , embedded in the hot gas (e.g. 
iMcNamara fc Nulsenll2007l ). The presence of bubbles is gen- 
erally taken to be a signature of AGN feedback. In this 
model, a fraction of the material cooling from the gaseous 
atmosphere is accreted by a supermassive black hole located 
in the central galaxy. This releases vast amounts of energy 
often in the form of outflows which couple to the hot atm o- 
sphere (e.g. IChurazov et al.|[20o3 : iBenson fc Babulll2009l ). 
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AGN feedback is widely considered to have important 
consequences for the evolution of single galaxies, as well 
as galaxy groups and clusters. For example, feedback is 
thought to be key in dete rmining the upper mass cut-off of 
the galax y mass function l|Benson et al.|[2003l : ICroton et al.l 
2006; Bo wer et al]|2006l ). balancing the radiative losses in 
cool-core galaxy clusters, as well as providing the non- 
gravitational 'pre-heati ng' that may be important for non 
cool- c ore c lusters (e.g. iBabul et al.l |2002| : [McCarthy et al.l 
|2004 [ 20081 ). Feedback also seems to play a role in deter- 
mining the relationship between the temperature and X-ray 
luminosity of hot a t mospheres across a ran g e of halo masses 



jBabul et al] l200d : iMcCarthv et al 
l2008l : iBower et al.l I2OO8I : iDave et al 



2004 iPuchwein et al.1 
2008 : |Popell2009t ).~As 



a direct consequence, feedback also regulates supermassive 
black hole growth (e.g. iSilk fc Reesi ilQQS : .Churazov et al] 
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l2005h and, therefore, its relation to the properties of the 
host galaxy. 

Considerable effort has been invested in understanding 
the impact of supermassive black holes. One of the key chal- 
lenges is understanding how the energy from the AGN af- 
fects and couples to the broader environment. Most studies 
that have sought to describe this relationship have focused 
on bulk motion, shock waves and pressure/gravity waves in- 
duced by AGN outbursts and the subsequent dissipat ion of 
the associated energy f see iMcNamara fc Nulsenll2007l . for a 
review). In this article we focus on the transport of material 
out of the cluster centre by AGN-blown bubbles. As an ex- 
ample, the filaments of cool material observed behind AGN- 
blown bubbles (e . g. Conselice et al.l I2001I : ICrawford et al.l 
I2OO5I : iHatch" et_alj|200y) are a strong indication of bubble- 
induced mass transport within several 10s of kiloparsecs of 
the cluster centre. The most obvious example is the Perseus 
cluster in which the filaments contain some 10* solar masses 
of cool material (e.g. Salome_et al. 2006, 2008). More gen- 
erally, some of the transport mechanisms we discuss in this 
paper may carry matter out to ~100 kiloparsecs from the 
cluster centre. 

Mass transport by bubbles not only reduces the amount 
of material available for forming new stars in the central 
galaxy, allowing the normally adopted stringent requirement 
that AGN heating perfectly balance cooling to be some re- 
laxed, but also affects the bubble dynamics and energetics. 
Hence, a better understanding of the main mass transport 
mechanisms is essential for describing the energy balance in 
the ICM. Here, we focus on three important processes by 
which bubbles can transport material: drift, wake transport 
and the entrainment of ambient material into the bubble. 

The aim of this article is to quantitatively describe these 
three main mechanisms and also the corresponding implica- 
tions. The discussion is largely analytical, intended to facili- 
tate a better understanding of the underappreciated aspects 
of these mechanisms and aid in the interpretation of both 
observations and numerical simulations. 

The article is arranged in the following way. Section 2 
outlines the basic bubble model that serves as a backdrop 
for subsequent discussions. Sections 3, 4 and 5 focus on each 
of the three mechanisms — drift, wake transport and en- 
trainment. In each section, we discuss a process as well as 
associated, potentially observable consequences, like trailing 
optical filaments. In Section 6, we consider how mass trans- 
port modifies the ongoing debate about whether or not AGN 
heating need balance cooling precisely. We summarize our 
key findings in Section 7. 



the end of the infiation process, the bubble is expected to be 
in pressure equilibrium with its surroundings. As a result of 
the mass density in the bubble being lower than that of the 
ambient medium, the bubble is buoyant and will rise, moving 
away from the AGN and in the case of galaxy clusters, the 
cluster centre, towards more distant regions of lower density 
and pressure. The velocity of the bubble is determined by 
the near balance between buoyancy and drag forces with the 
ICM. In response to dropping ambient pressure, the bubble 
will expand. The end stage of bubble evolution is not well 
understood. One possibility is that the bubble will be de- 
stroyed during the rise by hydrodynamical (Raleigh- Taylor 
and Kelvin-Helmholtz) instabilities; however, it is equally 
plausible that ambient conditions (magnetic fields on the 
bubble surface or viscosity in the ICM) suppress these in- 
stabilities, allowing the bubble to rise to a radius where the 
bubble density becomes equal to that of its environment, 
and buoyancy force vanishes. In the latter case, the bubble 
of relativistic fluid will hover at its equilibrium height as a 
"ghost cavity" , invisible in the X-rays though it may still be 
detectable in low frequency radio observations. 

In the context of mass transport, the most important 
stage in the evolution of the bubble is the second (or the rise) 
phase. The dynamics and thermodyna mics of this particular 
phase have beeri discu ss ed widely (e.g. Churazoy et al.ll2002l: 
[Enfilin fc Heinj |2002|: iMathews et all l2003l : iNusser et'al] 
[2006; Nul sen et allbOOTl T Typically, most analyses assume 
that the rise velocity of the bubble is sufficiently slow such 
that the response of the medium internal to the bubble to 
the bubble expansion is adiabatic. This approximation yields 
reasonable theoretical expectations. 

In addition, most treatments also assume that the bub- 
ble density is much less than that of the surrounding ICM 
and derive the relevant relationships in the limit of an ex- 
treme density contrast. For many of the calculations, this 
approximation is probably adequate. However, the density 
contrast determines the maximum height to which the bub- 
ble rises, and so is critical for determining the magnitude and 
the spatial scale of the mass transport phenomenon, as well 
as the efficiency with which energy can be extracted from 
the bubble. As we will discuss in the following sections, mass 
transport processes endow a bubble with an effective — and 
in the case of entrainment, a real — finite density contrast. A 
more complete description of bubble dynamics and thermo- 
dynamics must therefore allow for a finite density contrast. 
At relevant sections of this paper, we will therefore present 
modified relationships governing bubble dynamics and ther- 
modynamics that takes into account this aspect. 



2 BUBBLE MODEL 

The magnitude and spatial extent of mass transport associ- 
ated with buoyant bubbles depend heavily on the dynamics 
and the energetics of the bubble. It is, therefore, essential to 
have a model which accurately describes these processes. 

In broad brushstrokes, an AGN-inflated bubble is be- 
lieve to evolve as follows: an AGN jet inflates a bubble of 
what is generally assumed to be relativistic fluid in the intra- 
cluster medium. The details of the inflation process itself 
are not well understood and models for the bubble inflation 
range from nearly adiabatic to supersonic. In either case, at 



3 BUBBLE DRIFT AND TRAILING COOL 
FILAMENTS 

As a bubble rises, energy conservation requires that the 
change in its gravitational potential energy is accompanied 
by a gain in kinetic energy of the bubble and the ICM. The 
latter is induced by the ICM being pushed by the rising bub- 
ble, causing a net upward displacemen t behind the bubble. 
This phenomenon is known as "drift" (|Darwinll 19531 ) Q Fig- 

^ The phenomenon of drift is utilised in industrial mixing pro- 
cesses, called fluidiscd beds, which arc formed when gas is passed 
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ure 1 depi cts a schematic view of the bubble and its tra iling 
drift fe.g. ICrowdlioOsI : lYandliooal : ISchetz fc FuhsHigggl ). 

Drift is closely re lated to the bette r -known concept 
of "a dded mass" (e.g. iNusser et al. I I2OO6I : iPavlovski et al.l 
120071 '). a concept associated with moving bodies submerged 
in a fluid. The usual interpretation is that it represents 
the increased inertia associated with the work required to 
change the kinetic energy of the fluid flow around the body. 
The ambient fluid can be thought of as increasing the effec - 
tive mass of the moving object (e.g. lMilne-T hompsonl l 19961 ). 

In group and cluster environments, this trailing fluid 
can, under certain conditions, give rise to filaments of cool 
gas. In this section, we introduce the phenomenon of drift, 
discuss the associated characteristics, and explore whether 
the recently discovered optical filaments stretching between 
the AGN-inflated bubbles and the c luster centres (e.g. 
Conselice et al.l 120011: ICrawford et al.l l2005l : ISalome et al] 



20061 : iHatch et al.H2006l') are relate d to the drift. The kine- 
matic d^la~(e!gr^^cj^er|al] [200^) certainly suggests that 
the filaments may well be the result of cool material that 
was originally lifted from the cluster center and that some 
of this material is now falling back. 



3.1 Drift: a conceptual outline 

iDarwinI (|l953l ) considered the problem of a body moving 
uniformly through an infinite, inviscid, incompressible fluid. 
Importantly, Darwin studied the motion in the fluid frame of 
reference rather than the rest frame of the body in which the 
fluid particles follow simple streamlines. In the fluid frame, 
the trajectories are extremely complex, consisting of a large- 
scale looping motion and a permanent displacement in the 
direction of body's motion. 

The forward displacement occurs because a moving 
body pushes the fluid ahead of it causing the fluid to move in 
the same direction as the body. In due course, the forward- 
displac ed fluid s l ows d own, though in the situation consid- 
ered bv lDarwinI (|l953h never formally comes to rest behind 
the body. As an illustration, one can imagine a rising bub- 
ble approaching and traversing a stationary marked plane 
within the fluid column, oriented at right angles to its di- 
rection of motion. As the bubble passes through the plane, 
the fluid will be displaced upwards in the horn-like shape 
depicted in the third panel figure 1, with the fluid closer to 
the body experiencing a greater net translation. 

Strictly speaking, the fluid displacement has two contri- 
butions: Xd, a localised (positive) drift contribution which 
decays rapidly from the body, and Xr, a non-local (negative) 
reflux contribution which depends on the far-field boundary 
conditions and which decays slowly from the body. In the 
limit that the body starts infinitely far from the marked 
plane and then moves through the plane to infinity, the re- 
flux contribution is zero, so that the drift volume is single- 
valued. 

In the absence of gravity, the fluid displacement is per- 
manent and the length of the drift grows continuously while 



into a vessel containing a bed of solid particles. At a given gas 
injection rate, a stream of bubbles appears and rises through the 
bed with ta ils of material behind th em and thereby intensifying 
mixing (e.g. ICrowelbOOSi : lYanell2003l) . 



the total volume remains approximately constant. However, 
for drift behind buoyant bubbles in the intracluster envi- 
ronment, the stratification of the fluid cannot be neglected. 
Not only will gravity bring the drifting fluid particles to rest 
sooner, it can also potentially cause the fluid to fall back to 
its original position. As a result, in a stratified fluid, the up- 
ward displacement of material due to drift is always finite. 
The implications of drift in a stratified fluid are described in 
detail in the next subsection. Presently, we will assume, for 
simplicity's sake, that the effects of gravity can be neglected. 

The volume of fluid enclosed by the intial plane of 
marked fluid and the distorted plane is called the drift 
volume, Vdrift, and is depicted in figure 1. In an un- 
bounded flow, it is calculated from the integral of the per- 
manent displacement of fluid particles in the direction of the 
body's motion over the cr oss-sectional area of the flow (e.g. 
lEames fc Gilbertsonllioosi ). 

The displacement is given by the time integral of the 
relative velocity between the fluid in the drift and the 
fluid at infinity. It is calculated from the velocity poten- 
tial for a uniform flow of an incompressible , non- rotational 
fluid past a sphere, e . g. Landau fc Lifschitz j[l995l ) (see also 
IPavlovski et alj|2007 ': Dursi fc Pfr ommeil200a for recent ex- 
amples). Using this, [Darwin C1953i) showed that the drift and 
body volumes are related by a constant of proportionality, 
k, which depends on the geometry of the body; for a solid 
sphere, k — 0.5. In the case where b ubbles deform com- 
pletely into vortex rings, iDabiril 12006") found k = 0.72, so 
the value of k for an ordinary bubble should be somewhere 
between these limiting cases. For generality, we will use (e.g. 
lMilne-Thompsonlll996l : ISchetz fc Fuh^ligggl ) 



Vdrift = k{Vbub + Kv, 



(1) 



where Vwakc is the wake volume. (Fo r a full list of pa rameters 
and constants, see the Appendixra'l lDarwinI l|l953l ') was also 
able to show, quite generally, that the mass within the drift 
volume is equal to the added mass, so that 



Mdrift = fcpiCM(Vbub + Vwakc) ~ kM^in 



(2) 



where Mdis is the ICM mass displaced by the bubble. 

At the onset of the rise phase, the (initial) displaced 
mass can be estimated as Afdis.o = PoVq. Given the enthalpy 
of the bubble at the end of inflation and the beginning of 
the rise phase 

-Ebub.O = ' bubo _|_ PoVhubfi ~ — T -Po Vbub,0 , (3) 

ib — 1 tb — i 

where Fb is the adiabatic index of the bubble material, we 
can write 



M (Fb - 1) p.mp ^ms/' 
MdisA) = -; — —-Hbubfi ~ tU 



knTo 



( To 



1059 erg J\WK 



Mr. 



(4) 

Here, To is the ambient temperature where the bubble is 
inflated and we have assumed that the bubble contents are 
relativistic, hence Fb = 4/3. We note that the actual value 
of Fb will depend on the initial composition of the AGN 
outflow and whether the jet and the bubble inflation pro- 
cess entrains a significant amount of material from the ICM 
into the bubble. We will comment on this further when we 
consider entrainment. 

To estimate the drift mass that we would expect to be 
associated with an AGN-inflated bubble, we assume that 
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Table 1. Typical values of plicnomenological constants 



Parameter 


Plausible Range 


Best guess 


k 


0-1 


- 0.5 


q 


0-1 


~ 0.24 


Co 


> 0.45 


> 0.45 


a 


0-1 


< 0.05 



k = 0.5 since the bubbles tend to generally be more or 
less spherical. Accordingly, we would expect observations 
to show a total mass of a few x 10^ to ~ 10** M0 of material 
trailing behind typical bubbles. This value agrees reason- 
ably well with estimates of the mass in the recently discov- 
ered cool filament s that stretch betwe en the bubbles and the 
cluster centre fe.g lSalome et ah II2OO8'). 

In principle, it is possible to estimate k for an AGN- 
blown bubble directly from numerical simulations. Illustra- 
tive examples of the characteristic drift signature in galaxy 
cluster simulat i ons c an be found in Revaz et al ., (2008); 
iRoediger etHI l|200'H ), though we note that this effect can 
also be induced by circular vortex lines, which stretch 
bubbles into to ri (e.g. [Pavlovski et al■ll2007^ ■ Furthermore, 
IChurazov et al.l (I2OOII ) also estimated the mass of mate- 
rial uplifted by their simulated hydrodynamic bubble to be 
~ 10* Mq . This agreement is encouraging; however, we do 
caution that until we achieve a better understanding of the 
ICM microphysics (including a realistic treatment of mag- 
netic fields, thermal conduction and viscosity) any estimate 
for k based on numerical simulations is unlikely to signif- 
icantly improve upon the values quoted above. The same 
reasoning applies to three other phenomenological parame- 
ters used throughout this article: (1) q, the fraction of dis- 
placed mass trapped in the bubble wake; (2) Co, the bubble 
drag coefficient; (3) a, the entrainment coefficient. Plausible 
ranges and best estimates are given in table [1] 

Finally, we note that drift is enhanced at low Reynolds 
numbers, as can be se en in the viscou s hydr odynamic sim- 
ulations presented bv iRevnolds et al] (|2004l ). However the 
theory is much more complex as the drift volume does not 
have a well-defined maximum as well as being de pendent on 
boundaries in the system (see lEames et al]|2003l ). 



3.2 Observed cool optical filaments: a by-product 
of the drift? 

Optical filaments behind bubbles tend to be linear trails 
of cool material that extend for tens of kiloparsccs towards 
the c l uster centre (e.g. [C onsclicc ct al. 2001; Crawford ct al. 
I2OO5I : ISalome et all I2OO6I : iHatch et al.l 120061 '). These struc- 
tures are a strong indication of bubble-induced mass trans- 
port within several 10s of kiloparsecs of the cluster cen- 
tre. The mass associated with the filaments is estimated, 
in Perseus cl uster at least, to be ^ 10* solar masses of cool 
material (e.g. lSalome et al.ll2006l . 120081 '). 

The simplest explanation for the filaments is that AGN- 
blown bubbles lift cool material from the cluster centre when 
they detach and start to rise. The filaments would appear to 
grow when some of this material decouples from the bubble- 
drift system and falls back inward. Filamentary material 
close to the bubble would, in this picture, still be travel- 



ling upward with respect to the cluster centre, while ma- 
terial further away would be travelling towards the cluster 
centre. The filament would the n appear to be str etching, as 
indicated by observations (e.g. iHatch et al]|2006l ). 

The observable properties of filaments correspond quite 
closely to the general characteristics of drift. We, there- 
fore, investigate the circumstances under which these pro- 
cesses may produce flows that might be identifiable as fila- 
ments. Further to this, we also investigate additional pro- 
cesses which may explain the dumpiness of filaments on 
scales of tenths of kiloparsecs. 

As previously mentioned, the stratification of the ICM 
exerts a non- negligible infiuence on motions generated 
within it. The numerous different fiow regimes can be de- 
scribed by using the characteristic length scales of the sys- 
tem. Pertaining to drift behind buoyant bubbles, two useful 
scales are the Ozmidov and buoyancy lengths. The first pa- 
rameter relates to the maximum size of turbulent eddies in 
a stratified fiuid. Motions on larger spatial scales, therefore, 
cannot be turbulent alone and material will also be trans- 
ported vertically in the fiuid. Consequently, the Ozmidov 
length marks a characteristic scale of fiuid motion above 
which a filament may form behind a buoyant bubble. 

The scale length of vertical fiuid motion is given by the 
buoyancy length, which represents the complete conversion 
of vertical kinetic energy to potential energy in a stratified 
fiuid. It can be seen that filaments will be formed if the 
buoyancy length exceeds the Ozmidov length. These scales 
are considered in more detail below. 

To a reasonable approximation, the Ozmidov scale can 
be calculated in the following manner: the kinetic energy 
per unit mass required to displace a stably stratified fiuid 
a vertical distance Az is of order Az^N^, where A'' is the 
Brunt- Vaisala (buoyancy) frequency, while the energy asso- 
ciated with turbulence on the same spatial scale is (eAz)'^^'^, 
with e being the energy dissipation rate per unit mass. 
Equating the kinetic ene rgies gives the Ozmidov scale (e.g. 
iKantha fc Clavsonllioool ) 



where 



<; diner 
z d In z ' 



(5) 



(6) 



with g being the gravitational acceleration and a — PpiQ^ 
being the entropy index of the ICM. It follows that, in a 
stratified fiuid, there is a maximum kinetic energy associ- 
ated with the eddies that occur as a consequence of a given 
body's motion. Therefore, for objects that generate vertical 
disturbances on scales larger than Lq, the turbulent eddies 
cannot dissipate the kinetic energy. Instead, the majority of 
the kinetic energy goes towards increasing the potential en- 
ergy of the ICM by carrying denser material from closer to 
the cluster centre out to larger distances. Motion on much 
grea ter scales than Lq, is expected to generate gravity waves 
(c.f. lOmma et alll2004l l. In contrast, for fiuid disturbances 
that are smaller than Lq, turbulent eddies are highly effi- 
cient at disrupting the column of material displaced by drift, 
thereby preventing an appreciable change in the potential 
energy of the fluid - in this case no filament will form. 

The Ozmidov scale, therefore, marks the transition from 
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Layer of 



Wake- 




Figure 1. Schematic showing wake transport and drift due to passage of bubble. In panel (A) the bubble approaches an initially planar 
dye surface from below. (B) Distortion of the dye surface occurs as the bubble passes through the plane. (C) The volume between initial 
plane and horn-like distorted surface is the drift volume. 



motions that primarily generate turbulence to those that 
transport mass vertically. Since turbulence is associated with 
heating, the transition is equivalent to defining a threshold 
above which buoyant bubbles predominantly transport mass 
and below which they generate turbulence, which leads to 
heating. At this point it is also important to note that if 
the Ozmidov scale becomes comparable to the size of th e 
smallest eddies, the Kolmogorov scale (e.g. IChoudhuril 1998h . 
then turbulence cannot exist in the fluid. This is because the 
eddy kinetic energy is insufficient to circulate in the direction 
of the fluid stratification. However, unless the ICM viscosity 
exceeds the Spitzer value this is unlikely to be the case in 
clusters of galaxies. 

To simplify the expression for Lq, the turbulent dissi- 
pation rate for the bubble is taken to be e = Cturb'w'^ /r (e.g. 
IPennis fc Cha ndran 200^), where Cturb is a numerical con- 
stant, w is the bubble velocity, and r is the bubble radius. 
In order to describe the scaling properties of the Ozmidov 
length we will assume that the bubble ascends at its termi- 
nal velocity governed by the balance between buoyancy and 
drag, with a magnitude given by (e.g. Churazov et al.,,2001 ') 



^terminal 



(7) 



f8rg_Y 
\SCuJ ' 

where Cd is the drag coefficient. For brevity, equation ((7| is 
written tutorminai = (cterm^'S') aud the Ozmidov scale can 
be expressed as 



nurb^tcrm' 



din a 
din 2 



(8) 



where z is the bubble's displacement from the cluster centre. 

Equation (jS]) must be compared to the buoyancy length 
to determine whether filaments actually form. The buoy- 
ancy length is defined as the displacement of a fiuid parcel 
at which all its vertical kinetic energy is converted to po- 
tential energy. As before, the kinetic energy per unit mass. 



w , will displace fluid in a stably stratified atmosphere by 
a vertical distance Lbuoy , is of order L^^^^N'^ , so that (e.g. 
iKantha fc Clavsonlboodl ) the buoyancy length is 



w 
N' 



(9) 



For a bubble t ravelling at its term inal velocity w 
(ctcrmr-p)^/^ fe.g. IChurazov et aLlbOOll ). and 



-'buoy 



i 1 /dlna 



dlnz 



(10) 



Following the discussion in iDennis fc ChandranI ('2005'), our 
constant Cturb ~ 0.42/2 = 0.21. Then, using the fact that, 
for a spherical bubble Ctcrm = 8/ (3Cd), it can be shown that 
filaments can form behind buoyant bubbles provided 



-^buo> 
1^ 



1.7C^ 



I / r \ * / dlna\ * 



dlnz 



> 1. 



(11) 



Whether the criterion above is satisfied depends on the val- 
ues given to the phenomenological constants, Ctcrm and Cd, 
as well as observable parameters. It is particularly difficult 
to assign a characteristic value for the drag coefficient, Cd, 
since it depends o n both the Mach and Reynolds numbers of 
the flow (e.g. see lPope et al.l [20081 . and references therein), 
which are dependent on the nature of the ICM. Estimat- 
ing such parameters by comparison with numerical hydro- 
dynamic simulations is not necessarily representative. How- 
ever, a suitable guess can be obtained in the limit that Cd 
behaves as a conventional drag coefficient so that it varies 
from Cd ~ 30 for a Reynolds number of Re « 1, to Cd ~ 0.5 
for Re > 1000 (e.g. ICrowe et al.1 119981 ). This suggests that 
filaments are more easily generated in fiows with a low ef- 
fective Reynolds number, which may be the case for the 
rarefied ICM. Despite these uncertainties, equation pip re- 
mains informative because it indicates the functional rela- 
tionship between physical parameters. 
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The observational quantities are slightly better con- 
strained. In classical cool-core clusters, the entropy profile 
scales roughly li nearly with distance al l the way to the clus- 
ter centre (e.g. ISanderson et al.|l2009l ). with a logarithmic 
gradient close to unity. We also assume the empirical rela- 
tion that the bubble radius is linearly proportional to its dis- 
plac ement from the cluster centre, such that r/z ^ 0.3 — 0.5 
(e.g.lDiehletal ] |2008l : iMcNamara fc Nulsenll2007h . Combin- 
ing these values gives L^uoy/Lo ~ 1 — 3, indicating that 
material can potentially be displaced upwardly behind the 
bubble. Some of this material will come from the centre of 
the cluster and will be cool, meaning that an optical fil- 
ament could form beh ind the rising bubble. Interestingly, 
I Sanderson et al.l (|2009l ) present entropy profiles for a sam- 
ple clusters, of which at least A1795, A2199 and A478 are 
nearly linear entropy profi les and also show evidence for fil- 
amentary structures (see IPope et al.l |2008| . and references 
therein). We also note that the typical mass in observed fil- 
aments of ~ 10* Mq (e.g iSalome e t al. 2008) agrees well with 
drift mass estimates for typical ICM bubbles. Additionally, 
the closeness of I/buoy /^o to unity in clusters with nearly 
linear entropy profiles also suggests that buoyant bubbles in 
these systems do not generate gravity waves, though power- 
ful long-lived episo des of jet activity might (see, for example 
lOmma et alll2004l '). 

Classical cool core clusters, however, comprise only a 
small fraction of the clusters. In general, cluster entropy 
profiles tend to flatten close to the centre and then increase 
roughly linearly with distance from the cluster centre (e.g. 
iDonahue et al.l [20061 ). If the entropy profile has a flat core 
when the bubble begins to rise (as opposed to being ini- 
tially nearly linear and subsequently evolving into a cored 
profile due to the resulting heat injection), the logarithmic 
entropy gradient is approximately zero close to the cluster 
centre (see equation [8| and then roughly c onstant with a 
value near unity furthe r out, e.g. 1.0-1.3 (e.g. lDonahue et al.l 
I2OO6I : [Cavagnolo et al. I [2009). In this case, Lbuoy/io will be 
small near the cluster centre and the bubble motion will pri- 
marily generate turbulence which heats the central regions. 
Outside the central core, Lbuoy/io will take a value simi- 
lar to that quoted above for the classical cool-core clusters 
and some material may be displaced upwards by the bub- 
ble. However, since the upward displacement can only occur 
outside the central entropy core, there is unlikely to be a sig- 
nificant amount of cool material in the drift region. In this 
case, the drift will likely not appear as optical filaments. 



4 WAKE TRANSPORT 

Apart from the drift, a rising bubble can also transport 
material via "wake transport". When a deformable bubble 
starts to rise, it takes on a shape that is roughly spherical, 
but with an indentation at the bottom. This cavity between 
a true sphere and the actual bubb le shape with an indented 
base is known as th e wake (e.g. ICrowd [ioOSi : lYand [ioo^ : 
ISchetz fc Fuh3ll999l ). Figure 1 depicts a schematic view of 
the bubble, wake and drift. The "wake" in the present con- 
text is very different from the the turbulent wake that de- 
velops behind objects at high Reynolds numbers. 

Studies of bubbles and mass transport by bubbles in 
non-cosmological settings indicate that line vortices gener- 



ated behind the bubble as it begin to rise lift, trap the am- 
bient medium in the vicinity (i.e. material in the cluste r 
centre) within the "wake" cavity (c.f. iPavlovski et al.ll2007l '). 
This trapped material is then carried up by the bubble. 

In the first instance, material in the wake increases the 
average density of the bubble and as a direct consequence 
reduces its ascent velocity and the maximum height to which 
the bubble rises. In addition, trapped material can also drop 
out of the wake, and this material will further contribute to 
the formation of filaments behind AGN-blown bubbles. 

4.1 Wake: a conceptual outline 

Studies of fluidised beds indicate that a typical wake con- 
tains a fraction q of the initial mass displaced by the bubble: 

Mwakc ~ gMdis.O, (12) 

where q « 0.24 and Mdis.o is given by equation (Q. As- 
suming that this laboratory value of q is applicable to the 
intra-cluster medium, the above relationship suggests that 
AGN- inflated can easily Uft few x 10^ to ~ 10* Mq in their 
wakes. 

In addition to facilitating the transport of material out 
of the cluster centre, the wake — and specifically, the mass 
in the wake — increases the effective average density of the 
bubble (though it does not affect the bubble's actual density 
contrast as measured from X-rays). As we discuss below, this 
modification impacts not only the maximum height that the 
bubble can ascend to, but also the efficiency with which 
energy is extracted from it. 

We do, however, note that as in the case of the geometric 
constant k in the estimate for the drift mass, we acknowledge 
an uncertainty: it is unclear whether the laboratory value of 
q is applicable to the intra-cluster medium. It is difficult to 
propose a better motivated value because our understand- 
ing of the ICM microphysics is highly incomplete. This also 
means that the estimates of q based on numerical simula- 
tions are not necessarily more believable than the laboratory 
value deduced from the observations of fiuidised beds. On 
the other hand, one would expect that q does not exceed 
unity, since this would imply that the wake contains more 
mass than originally displaced by the bubble; consequently, 
unless q 1 for the ICM, our estimates for the wake mass 
in the range of ^ 10* M0 should be reasonable. 

4.2 The implications of wake transport 

Studies of bubbles indicate that dur ing the ascen t, the bub- 
ble wake can lose blobs of material (|Crowdl2005l ). Formally, 
this happens if the terminal infall velocity of a blob exceeds 
the local upward velocity behind the bubble. Once detached 
from the wake, the subsequent evolution of the material is 
unclear. In the simplest picture, the blob simple falls back 
inward. The presence of the ICM complicates the scenario in 
that the associated ram pressure will likely strip and stretch 
the blob into filamentary structures. Moreover, the material 
from the shredded blobs may also get caught in the drift 
behind the bubble and some of this material may s till get 
displaced to larg er distances from the cluster centre (| Crowd 
l2005l : lYang||200^ ). 

This discussion would suggest that both drift and leak- 
age of material from the bubble wake can contribute to the 
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formation of trailing filaments. However, the bulk of the 
wake material will tend to fall towards the cluster centre. 
As a result, the trailing material is unlikely to exhibit a sin- 
gle characteristic kinematic signature. Rather, the kinematic 
structure will depend on the relative importance of drift and 
leakage from the wake, and will probably vary significantly 
over time as different processes become prevalent. This issue 
bears further investigation. 

As already noted, both the wake and drift affect the 
dynamics of the bubble, but in slightly different ways. To 
illustrate this, the equation of motion of a bubble carrying 
a wake, and subject to the action of drift, can be written as 

d[(Mbub + Mwake + Mdrift)w] 



dt 



(13) 



A/dis 



where Mdis is the mass displaced by the bubble, which falls 
with time, and F is the drag force 

„ 1 



-CoApicMW 



(14) 



and A = Tvr is the cross-sectional area of the spherical 
bubble. 

From the left hand side of equation psp it can be seen 
that a bubble of constant mass, Mbub, which transports the 
additional masses, Mwake + Afdrift, accelerates as if it has a 
larger effective mass of Mbub + Mwake + Mdrift- However, in 
the limit that the acceleration is zero, the bubble behaves 
slightly differently again; the terminal velocity of the bubble 
described by equation p3p . can be expressed as 



tt^tcrminal — 



8gr_ 
3Cd 



A'/w 



+ Mb 



Md 



(15) 



Equation (|15|) can be interpreted as the terminal velocity of 
a bubble with an effective density contrast of 

Afdis 

^ M Tm • ^^^^ 

Mwake + JWbub 

Denoting the effective density contrast when the bubble 
first starts to rise after inflation as ?7cft,o, we note that 

Vo 

VcS,0 ~ 



where rjo = po/Pbub,o = A/dis, o/Mbub is the actual initial 
density contrast and Mdis.o = PoVo is the initial displaced 
mass. Here, po is the ICM density where the bubble is in- 
flated, and Pbub.o, Vo are the bubble mass density and the 
bubble volume once the bubble is inflated. A typical value of 
the actual initial density contrast taken from the literature 
is ?7o ~ 100 fe.g. iGhizzardi et al.ir200 4''). With mass in the 
wake characterized by q « 0.24 (see [Vang 2003), r/cff.o ~ 4. 
In this derivation, we have assumed that the mass of the 
bubble remains constant during the rise. This is equivalent 
to assuming that regardless of whether the ICM is entrained 
within the bubble at formation, the bubble does will not sub- 
sequently entrain any additional material during the rise. We 
have also assumed that the wake mass remains constant. The 
displaced mass is allowed to change since both the bubble 
volume and the ambient density vary with distance from the 
cluster centre and for completeness, we note that this mass 
will decrease as the bubble rises. 

According to equation (|15p . the low effective density 
contrast reduces the terminal velocity of the bubble. In the 



initial stages, r/ofi.o ~ 4, results in a lower terminal velocity 
by a factor of ~ 0.9. And even though the effect is not large, 
it is worth bearing in mind since the bubble rise timescale, 
a quantity commonly used to estimate AGN power, is in- 
versely proportionate to the velocity. 

A lower density contrast also reduces the equilibrium 
height to which the bubble can rise and a sufficiently low 
contrast can result in the bubble not being able to rise be- 
yond the very central regions of the cluster. This can be 
illustrated as follows. 

Under the simplest circumstances, and ignoring fluid in- 
stabilities, the buoyant bubble will rise through the stratified 
atmosphere and expand adiabatically to maintain pressure 
equilibrium with its surroundings. Consequently, the bubble 
density can be r elated to the ambie nt density and tempera- 
ture by (see also lHinton et al"]|2007l ) 



f PiCM 
Pbub,0 V Po 



Pbub 



Tk 



To 



(18) 



where To and po are the ICM temperature and density, re- 
spectively, where the bubble was inflated, Ticm and picM 
are the ambient temperature and density at the bubble's 
current location. 

This relationship shows that the bubble density drops 
more slowly than the ambient density and there will be a 
location where the densities become comparable. In the ab- 
sence of a wake, the buoyancy force would also vanish at this 
height, and the upward motion of the bubble would stagnate 
since this motion is driven by buoyancy. In practice, due to 
the kinetic energy gained by the bubble, it will not immedi- 
ately come to halt at the equilibrium height. Instead, it will 
overshoot this location. Above the equilibrium height, the 
bubble will be denser than its surroundings and experience 
a gravitational restoring force towards the cluster centre. 
Consequently, it will oscillate around the equilibrium radius 
at a pproximately the Brunt- Vaisala (buoyancy) frequency 
(e.g. IChoudhurilll998l ). The oscillation will, in due course, 
be damped by the drag force but while it is oscillating, the 
bubble will generate fluid disturbances that travel towards 
and away from the cluster centre. Consequently, a fraction of 
the kinetic energy gained during the ascent can go towards 
heating the central regions of the ICM. 

The arresting criterion for the bubble- wake system, in 
the presence of the wake, is obtained from equation (|13p by 
noting that at the equilibrium height, Zcq 



Mw 



^(Zeq) ~ A/dis (^cq) " Mbub 



(19) 



where Mdis(zeq) = piCM{zcq)Vbubizc(i) is the displaced mass 
at the equilibrium location. This criterion is equivalent to 

'7eff(Zeq) = 1. 

Since T^off can be written as 



VoS = ??cff ,0 



f piCM \ / Tlc^ 



\ po J V Fo 
then 77cft (zoq) = 1 at the equilibrium height implies. 



pICM(-Zeq) 

po 



■ VeS,0 



7lCM(2cq) 



To 



(20) 



(21) 



where TicM(zcq) is the ICM temperature at the equilib- 
rium height. For coo l-core clusters w e typically expect 
ricM(^eq)/ro ~ 3 (e.g. lAUen et Ikl bOOll ). 
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Modelling the ICM density distribution using a /3-profile 

2 1-/3 



picmiz) ~ PiCM.o 



1 + 



(22) 



V 2ICM,0 , 

where picM,o is the central ICM density, zicm.o is the ICM 
scale height, and z is the distance from the cluster centre, we 
can estimate the equilibrium height to which an adiabatic 
buoyant bubble carrying a wake, inflated at height zq, will 
rise in such an atmosphere: 

(23) 



1 + 



2lCM,0 



1) 




^<(^b-l) 


f ZlCM,0^ 




[ To \ 







^ VcS,0 



Restricting ourselves to classical cool core clusters that 
generically exhibit very small scale heights and hence, nearly 
power-law density proflles, we flnd that 



'/cff.O 



Ticm{zc 



To 



3(rb-i) 



(24) 



In Figure 2, we plot equation (|24|) as a function of g = 
Afwako/Mdis,o for three different values of the initial density 
contrast 77efi,o = 1000,100,10,4, assuming Fb = 4/3, P = 
3/4 and TicM{zcq) /To — 3 (as is appropriate in cool core 
clusters) . 

It can be seen that a low effective density contrast due 
to the wake can significantly affect the equilibrium height, 
resulting in the bubble being trapped closer to the cluster 
centre. Specifically, we find that if we ignore the wake (i.e. set 
q = 0), then T^cff.o = Vo ^ 100 and Zcq ~ few x 10'^ zq, which 
probably greatly exceeds the cluster virial radius even if zo 
is small (i.e. zo ~ 1 — 30kpc). For the same situation, if 



0.24 and r/o = 100, then ?7eff,o ~ 4 and 2e 



4.5zo. If 
= lOkpc, 



the bubble is infiated and begins its ascent from zo 
it will arrest at ~ 50 kpc. 

For completeness, we note that the ratio of the final to 
initial pressure is 



To 



(25) 



TlCM(Zcq) 

It can also be shown that the bubble's density contrast 
is important for governing the efficiency with which energy 
is extracted. The total energy available for heating the ICM 
and transporting mass is obtained by integrating the buoy- 
ancy force acting on the bubble along its path of motion (c.f. 
IChurazov et~al]|2002l : iNusser et ai]|2006l ') 



/ -Fbuoydz = / Vhuhpbubgdz 
J zn J zn 



(26) 



VhubPlCMgdz 



Afwakc^dz 



where g is the gravitational acceleration. Making use of the 
fact that the ICM is in hydrostatic equilibrium, the second 
term on the right hand side of equation (|26[) can be re- 
written as 

dP 

VbubPiCMgdz Vbub-j— dz = VbubdP, (27) 



where P is both the pressure in the bubble and the ICM 
pressure. These two are assumed to be equal since the bub- 
ble is in pressure equilibrium with its surroundings. Using 



equation (|25p and the fact that the bubble expands adiabat- 
ically, 

r-P(Zoq) p 

VbubdP = ,^ \, PoVo 



(rb - 1) 



7lCM(Zcq) 



'?cff,o2o 



1 



(28) 



where we have used equation (|25p to derive the latter rela- 
tionship. 

The remaining two terms on the right side of equation 
H26|) can also be approximated if we make some reasonable 
assumptions about the bubble-wake-ICM system. For ex- 
ample, for a bubble- wake system where Mbub + M„akc is 
2 constant, the requirement that the ICM is in hydrostatic 
equilibrium allows us to write 



VbubPbnbgdz+ I Afwakeffdz (29) 


Mdi.,0 /■^'"""^ dP 



?7cff,o 



PiCM 



To make further progress, we restrict ourselves to classical 
cool core clusters whose entropy profiles sca le nearly linearly 
with distance from the cluster centre (e.g. ISanderson et al.l 
liooi) and in keeping the observations, we approximate that 
the ICM density profile for classical cool core clusters as a 
power-law: picM ~ po {z / zo)~'^^ ■ This profile is equivalent 
to the /^-profile in equation \22\ in the limit zo S> zicm,o, 
where zo is the radius at which the bubble is inflated. For 
/? = 3/4, we find that 



Mwakc gdz 



Fb 



VbubPbub gdz+ 



where we have applied equation 

Combining equations (|28[l and (|30[l . we find that 



(30) 



PoVo ( Ticyi{za<i] 



(Fb — 1) r;off,o V ^cff.oTb 



f 

J Zn 



(rb - 1) 
1 



PoVb 

TlCM(Zeq) 



(31) 



In 



??cfl,0 



rjcRfiTo 



l(2eq) 



For ricM(zGq)/ro ~ 3, comparing the energy available for 
heating when r^oft.o ~ 100 and ryeff.o ~ 4, which is the more 
appropriate value when the wake mass is taken into account 
and corresponds to g ~ 0.24, we find that heating is re- 
duced by a factor of '--^ 6 in the latter case. Physically, this 
means that some of the bubble's energy goes into raising 
the potential energy of its content, including the wake, and 
is therefore unavailable for heating the ICM. The lower the 
density contrast between the bubble and its surroundings, 
the lower the fraction of energy released during the ascent. 



5 ENTRAINMENT 

Mass transport may also proceed by the entrainment of ICM 
material into the bubble — the drawing-in and incorporation 
of material from the ICM into the bubble body — and its 
subsequent displacement to larger radii in the cluster, as the 
bubble rises (c.f. iMcCarthv et al... 20081 

Entrainment can impact bubble dynamics in two impor- 
tant ways: (1) It causes the bubble to expand more rapidly 
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Figure 2. Plot showing the equilibrium height as a function oi q = M,vako/^dis,Oi for three values of rjo using equation I I24I I. For this 
plot, Tb = 4/3, 13 = 3/4 and TicM(zeq)/ro = 3. 



than predicted by a simple adiabatic model, and (2) it lim- 
its the height to which a bubble can rise. The importance of 
this particular channel, however, depends on the extent to 
which material is drawn into the bubble. 

Establishing the extent of entrainment in AGN-inflated 
bubbles is a non-trivial task. Entrainment is an extremely 
complex phenomenon and can occur in several ways. For ex- 
ample, the ambient material may not be displaced by the jet 
with 100% efficiency, resulting in the entrainment of some 
ICM material within the bubble. Kelvin-Helmholtz instabil- 
ities due to shear between the jet and the ICM can also draw 
ambient material in through the sides of the jet. Addition- 
ally, Rayleigh- Taylor instabilities on the upper bubble sur- 
face can allow the ICM material to penetrate, resulting in a 
time- varying mass entrainment rate as the bubble rises. On 
the other hand, given that the observed bubbles appear to be 
highly resilient to instabilities that can transform spherical 
bubbles into vortex rings, a number of researchers have sug- 
gested that the AGN-blown bubbles may either be draped 
by surface magnetic fiel ds, or that the ICM may be more vis - 
cous than thought fe.g. lDe Youn3[2003 : Kaiser et al.ll2005l ). 
Both surface magnetic fields and enhanced ICM viscosity 
dramatically inhibit mixing and the growth of fluid instabil- 
ities. But even in these situations, intra-cluster material can 
still enter the bubble. For example, in the case of a bubble 
surface draped by surface magnetic fields, the ICM can still 
leak into the bubble as a consequence of magnetic reconnec- 
tion events and diffusion across the bubble/ICM boundary, 
(see .Pope 2010). 

The upshot here is that entrainment of the ICM into 
AGN-blown bubbles is not governed by hydrodynamic pro- 
cesses alone. As a result, findings from current generation of 
numerical hydrodynamic (and even magnetohydrodynamic) 
simulations are not directly applicable; these do not, as of 
yet, adequately capture the necessary microphysical pro- 
cesses. This then argues in favour of treating the process 
in a more general manner, as we attempt to do here, with 



explicit (rather than implicit) assumptions about the nature 
of the fluid flow. 



5.1 Entrainment: a conceptual overview 

As a starting point, we turn to the conservation equations for 
mass, momentum and thermal energy that are often applied 
to the study of plumes in the terrestrial atmosphere (e.g. 
iTurneilflQea ). In doing so, we avoid having to make any im- 
plicit assumptions about the fluid behaviour of the ICM, in 
contrast to standard hydrodynamic simulations. For exam- 
ple, the flow of an inviscid, rarefled fluid around an obstacle 
can have a low effective Reynolds number (e.g. _Shen 2005i), 
which cannot be accounted for in the continuous fluid ap- 
proximation. Furthermore, the Reynolds number range ac- 
cessible to hydrodynamic simulations is Re < n^, where n 
is the number of computational cells across the width of 
the flow. Consequently, coarser spatial resolution implies a 
larger effective viscosity that may act to reduce the growth 
rate of instabilities on the bubble surface, which affects the 
entrainment rate. Instead, we model the effective Reynolds 
number of the flow, affected by rareflcation and viscosity, 
through the drag and entrainment coefflcients. In partic- 
ular, the dr ag coefflcient is related to the Reynolds num- 
ber (e.g. see iPope et al.ll2008l . for an approximate functional 
form) and hence the effective viscosity. Furthermore, we cal- 
ibrate the entrainment rate to match mixing rates due to 
the growth of Rayleigh- Taylor instabilities in hydrodynamic 
simulations of buoyant bubbles in cluster at mospheres (e.g. 
IChurazov et all l200ll : iPavlovski et all l2008f ). Applicability 
to the ICM is also extended by accounting for the different 
adiabatic indices of the bubble and the ICM. 

Integrating the equations then yields the bubble radius, 
velocity and density as it evolves in an external medium. 
Importantly, since the bubble is assumed to be spherically 
symmetric, the integration can be performed to very high 
precision so that numerical errors are significantly dimin- 
ished. It is also staightforward to ensure the equations are 
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entirely self-consistent and that the appropriate quantities 
are conserved. 

Nevertheless, certain assumptions are required to solve 
the conservation equations given below; the most impor- 
tant being that the bubble remains spherical, intact and 
in pressure equilibrium with its surroundings. This is ob- 
servationally justifiable; with the possible exce ption of M87 
which appears to exhibit a vortex ring structure l|Owen et al.l 
The ICM is also taken to be hydrostatic and we omit 
any effects due to radiative processes and thermal conduc- 
tivity. Other than this, the derivation does not make any 
assumptions about the shape of the ICM density or tem- 
perature profiles. These assumptions are similar to many nu- 
merical simulations of bubbles and are permissible for the 
following reasons. Radiative losses from the bubble depend 
on its particle content, which is difficult to determine and 
might plausibly be small. Thermal conduction of heat across 
the bubble/ICM must be strongly suppressed, probably by 
magnetic fields, other wise bubbles would rapidly evaporate 
l|Pavlovski et al.ll2008l ). 

To make the model as transparent as possible, we have 
employed the well-known entrainment hypothesis which 
states that the velocity, v, at which material joins the bub- 
ble is proport ional to the ascent velocity of the bubble, w 
(|Turneij|l969l ). The constant of proportionality, a, is called 
the entrainment coefficient and is defined by 



(32) 



In a uniform medium, and in the Boussinesq limit, the typi- 
cal entrainment coefficient for a momentum-driven o utflow is 
Q ~ . 05 and ~ . 08 for a buoyancy-driven flow (e.g. lTurneil 
ll986l ). lDe Yound (|2006[ ) also describes how the entrainment 
coefficient may vary with the density contrast, rj, and Mach 
number, Mo, such that a ~ 0.2(Ma/0.5)"^ (r;/2)^/^ How- 
ever, this functional form is representative of mixing due 
to non-linear Kelvin-Helmholtz instabilites, which are sup- 
pressed on the surface of a magnetised AGN-blown bubble. 
As a result, we assume a to be a constant. 

Following equation (|32p . the equation of mass conser- 
vation for an entraining bubble is 



dMb 



At 



= QflbubpiCMlU, 



(33) 



where Ob is the surface area of the bubble. The relation is 
best expressed in terms of the distance from the cluster cen- 
tre, so, assuming that ^ = f^g^ and the bubble is spherical, 
we can write 



dMbub , 2 

— ■ = inr apicM. 

az 



(34) 



Substituting A/bub = PbubVbub into equation (|34|) gives 



dr _ ^ picM r I dpbub 
Az pbub 3 pbub Az 



(35) 



Using equation H35p . forms the basis for calibrating a against 
numerical simulations of bubbles in cluster atmospheres. 
The length scale over which entrainment dominates the mass 
of the bubble can be defined as the distance, L, during which 
the bubble entrains a mass comparable to the initial dis- 
placed mass. To first order, the entrainment length is given 



by 



L : 



dMb 



Az 



11 



5 kpc 



a 
(U5 



kpc 



(36) 

where ro is the initial bubble radius. □ Hydrodynamic sim- 
ulations of bubbles in galaxy clusters show that bubbles al- 
most completely mix with their surroundings before they 
have ascended a distance o f roughly 2-3 init i al bu b ble radii 



(e.g. IChurazov et al.l 



I2OO7I : IPavlovski et al.l 



200fl : iRobinson et~all |2004| : iGardinil 



20081 ). On this evidence, the effective 



entrainment coefficient due to the development of Rayleigh- 
Taylor instabilities must be a? ~ 0.1 — 0.15. The effect of 
viscosity and magnetic fields are equivalent to a reduction 
of a. 

To proceed further, it is necessary to consider the mo- 
mentum and thermal energy equations. The momentum 
equation is given by equation (fT3)) using fepicMVbub = Mdrift 
and a constant wake mass of gAfdis.o- The thermal energy 
equation can be derived from the following statement: the 
addition of a mass of gas, dA/bub, at temperature Ticm, and 
at constant pressure, changes the internal energy of the bub- 
ble, Cp^bubTbubMbub- In this description, Cp,bub is the specific 
heat capacity at constant pressure, per unit mass of the 
bubble and Cp^bub = rbfcB/[p'Tip(rb — 1)], where Fb is the 
adiabatic index of the bubble material. The thermal energy 
equation is then 



CpJCMllCMdAfbub + Cp,bubrbubA'/bub = 
(Tbub + drbub)(A/bub + d A/bub )Cp, bub, 



(37) 



where Cp,icM = FicMfcB/[p»7ip(FicM — 1)], is the specific heat 
at constant pressure per unit mass of the ambient gas. 

Expanding the right hand side of equation (|37|) up to 
first order terms and rearranging gives 



dpb 



dA-fb 



Pbub 



Afb 



riCM(Fb — 1) Pbub 



rb(FlCM — 1) pICM 



(38) 



The density of the bubble is also affected by adiabatic expan- 
sion. Assuming the atmosphere is hydrostatic, the change in 
density can be written in terms of the gravitational acceler- 
ation of the cluster 



dpb 



Pbub 



AP 

IVP 



firripgAz 

FbfcBriCM 



(39) 



However, because the bubble now contains a mixture of fiu- 
ids with different adiabatic indices, its expansion must be 
described in terms of a single, time-dependent effective adi- 
abatic index calculated in the following way. Entraining a 
mass of material, AM, at the same pressure, is equivalent 
to entraining a volume AV^ = AM/picM- This increases the 
internal energy of the bubble by an amount PAV/ (Ticm — I) , 
after which the total must be 



P(Vbub + AV) 



Vb, 



AV 



(40) 



where F is the effective adiabatic index. Then, after the ith 



^ The critical value of a above which equation H36| l strictly holds 
occurs for L = zjcMfi, so that Ocrit = ''o/(32/CM,o)- For the 
current setup with an initial bubble size of ro = 5 kpc, and a 
density scale height of zjcMfi = 18.6 kpc, we find Ocrit ~ 0.09. 
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timestep of the numerical integration, the updated value of 
r must be 



Ti+i - 1 iVi + AVi) 



Vi 



AVi 



^Ic^ 



(41) 



where AVi is entrained volume during the ith timestep, so 
that Vi + AVi is the total bubble volume after the timestep, 
Fi was the effective adiabatic index at the beginning of the 
timestep and Ficm is the adiabatic index of the ICM. In 
the subsequent equations and numerical calculations we use 
Fb = Fi. 

Equations ((33, (O, (|3HI and can then be rear- 
ranged to yield simultaneous differential equations for the 
bubble radius, density and velocity 



dpb 



dz 



dr 
di 

3a 
r 



Ficm (Fb - 1) tJ^rriprg 

t — 

Fb (Ficm — 1) SFb^BTici 



PiCM 



and 

dw 
dl 



-3awpicM/r 

Ptot 



1 + 



(ricM-i)^'"'_ 

^ A"^pPbubg 
FbfcBTlCM 



FicM (Fb — 1) 



(42) 



(43) 



+ 



Fb (Ficm — 1) 

9 (ptot - (1 + fc)pICM) 

W ptot 

(Fb - l)fimpwg fcpicM 

Fb^BTlCM Ptot 



(44) 



VbubWptot 

where Vhub is the bubble volume, fcpicM Vbub is the drift mass 
and defining a total effective density 

qMdis.o 



Ptot = Pbub + kpiCM 



14, 



(45) 



with qMdisfi being the wake mass. 

The terms in equation (|44p represent entrainment, 
buoyancy, added mass and drag, respectively. A further brief 
inspection reveals that the adiabatic relations are recovered 
in the limit that entrainment is unimportant (a — >■ 0). The 
effective adiabatic index of the fluid mixture within the bub- 
ble is given by equation (|41[) . 

In addition to these equations, we must also spec- 
ify functional forms for the ICM density and temperature, 
which can be used to calculate the background gravitational 
field. As an example, the ICM density can be modelled as a 
/3-profile, e.g. equation (|22p . For t he temperature profile, we 
use the Gaussian function used bv lChizzardi et aU (|200J) to 
describe the ICM temperature of the Virgo cluster 



TiCM ^Ti ~T2 exp 



(46) 



where Ti is the ICM temperature at large radii, and Ti — T2 
is the temperature in the centre of the cluster, zt is the scale 
height of the temperature profile. 

The background gravitational acceleration of the cluster 
potential is determined from the temperature and density 
profiles by assuming the gas is in hydrostatic equilibrium 



kB ( dTlCM , TiCM dpicM 



/imp \ dz picM dz 



(47) 



W ithin this framework, we take a fit from IPope et al.l 
(|2006l ) to the density profile of t h e Hyd ra cluster based 
on the data given by iDavid et al.l (j200lh . The values are 
representative of general cluster density profiles and are 
po = 0.7 X 10"^^gcm-^ (3 = 0.72 and zicm.o = 18.6 kpc. 
The temperature is taken to be a typical cool core cluster in 
which the central temperature is a factor of 3 lower than at 
large radu: Ti = 3 x lO'' K, T2 = 2 x lO'' K with zt = 5.0 kpc. 

To illustrate the model, we have numerically integrated 
equations (|42ll . (|43|l and (|44|l using a standard fourth-order 
Runge-Kutta method. For the initial conditions, the bubbles 
were placed 10 kpc from the cluster centre, with a radius of 
5 kpc and a density contrast of ~ 40. In both entraining 
and non-entraining cases, the added mass coefficient was set 
to fc = 0.5, with q = 0.24. The initial adiabatic index of 
the bubble is Fb,o = 4/3, Ficm = 5/3 and Co = 0.5. By 
definition, a = in the non-entraining case, while for the 
entraining bubble we set a — 0.15, as expected for mixing 
due to Ray leigh- Taylor instabilities. 

Figure 3 shows the density contrast of bubbles evolving 
in the Hydra-like cluster atmosphere. In these plots, the solid 
line represents the non-entraining bubble, while the long and 
short dashed lines indicate the behaviour of the entraining 
bubble with and without a wake, respectively. 

As can be seen from figure 3, it is clear that entrain- 
ment increases the rate of expansion of the bubble over and 
above what would be expected due to adiabatic expansion. 
According to equation (|42|l . for large values of a the half- 
opening angle of a bubble would tend to 



tan(e) = ^ ~ ^a, 



(48) 



using Fb = 4/3 and Fk 



5/3. 



5.2 Entrainment: Observational implications 

As already discussed in the context of the bubble wake, 
bubble dynamics and energetics both depend on the den- 
sity contrast of a rising bubble. Entrainment obviously con- 
tributes to the lowering of this contrast. In addition to this, 
understanding the impact o f entrainment is of particular in- 
terest because iDiehl et al.l ((2OO81) have noted that bubble 
radii are approximately half their distance from the clus- 
ter centre. They argue that this correlation is too large and 
too constant to b e attr ibutable to simple adiabatic evolu- 
tion. iDiehl et al.l (|2008l ) favour a model where bubbles ar e 
inflated by cu r rent-d ominated jet (|Nakamura et al]|2006l ). 
iBrfiggen et al] l|2009l ), however, argue that this behaviour 
may well be explicable by a combination of projection effects 
and modifying the standard adiabatic model to include the 
influence of entrainment. 

According to observations, the typic al opening angle o f 
a bubble is roughly tan(6') ~ 0.5, (e.g. iDiehl et al.1 120081 ). 
If entrainment is the physical process which leads to this, 
equation (|48p implies that a ~ 0.8. This, in turn, implies 
an entrainment length-scale of L ~ ro/2. In other words, 
the bubble could only rise a distance equivalent to the ini- 
tial bubble radius before it became indistinguishable from 
its surroundings. It is difficult to rule out this possibility en- 
tirely since we cannot know the bubble radii at earlier times. 
However, bubbles clearly pers ist out to larg e dista nces - the 
cavity statistics provided by iBirzan et al.l l|2004l ') indicate 
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Figure 3. Plot showing the density contrast for the entraining bubble (a = 0.15,0.05) and a non-entraining bubble (o = 0). Note that 
the density contrast rapidly becomes small in the entraining cases. In fact, the final density contrast is less than unity because is 
always less than Ficm, see equation II43I I. 



that there are bubbles located more than 20 kpc from the 
cluster centre in Perseus, A4059, A133, A2199, RBS 797, 
A2597, Hydra A, Cygnus A. 

Moreover, ISanders fc FabianI (|2007l ) argue, based on an 
analysis of the bubbles in Perseus, that the volume fraction 
of thermal material in bubbles cannot exceed 50%. The con- 
version of the volume fraction into a mass fraction depends 
requires making certain assumptions about, for example, 
the temperature of the entrained material and whether it is 
heated significantly as it enters the bubble. Nonetheless, it is 
possible to derive some general expectations: if the entrained 
material is at the ambient temperature then the density con- 
trast of the rising bubbles must be high and the entrain- 
ment process must be very inefficient (i.e. oi must be small), 
otherwise the bubbles would not ap pear as depressions in 
the X -ray surface brightness mapsfe.g. lMcNamara fc NulsenI 
I2OO7I ). If the entrained gas is heated above 15 keV, it would 
be difficult to detect using existing X-ray telescopes and 
therefore, the appearance of AGN bubbles as depressions in 
the X-ray surface brightness maps offers no constraint. On 
the other hand, if the bubble and the ICM are in pressure 
equilibrium, then the fraction of the bubble mass attributed 
to entrainment depends roughly on the ratio of the ambi- 
ent temper ature to that of the new ly heated material in the 
bubble fe.g. lMcCarthv et al.ll200^ \ If the entrained material 
is heated significantly, the entrained mass can again not be 
very large. For example, if the entrained material is heated 
to be, say, five times the temperature of the ambient gas, the 
bubble will contain at most one fifth of the mass originally 
displaced by the inflation process. 

Consequently, the upper limit of 50% for the volume 
fraction of the thermal content strongly suggests that the 
bubbles cannot be entraining rapidly and the large opening 
angles must have a different explanation. This can be illus- 
trated using the similar reasoning used to derive equation 
l|40|l . Specifically, the pressure fraction, /i+i, contributed by 
entrained material after the ith timestep of the numerical 



integration using 



/i 



+1 



Tb.o 



1' 



(49) 



where Fb.o is the initial adiabatic index of the bubble be- 
fore it entrained any material from its surroundings and 
/i+i = thermal pressure/total pressure. Note that this is the 
volume-averaged pressure fraction. Therefore, since the pres- 
sures of each component must have been equal, and given 
that material is entrained at constant pressure, the volume- 
averaged pressure fraction is equivalent to the volume frac- 
tion of thermal material. 

Figure 4 clearly shows that, for a — 0.15, the thermal 
pressure fraction grows rapidly to greater than 0.6 before 
the bubble has risen 10 kpc from its original location. How- 
ever, for Q = 0.05, the pressure fraction is 0.5 at roughly 
40 kpc from the cluster centre and rises much more slowly. 
For a = 0.01, the pressure fraction remains <50% even be- 
yond 150 kpc. Consequently, a thermal volume fraction of 
~ 50% suggests either that the thermal content was injected 
(or entrained) during the initial inflation, or that the sub- 
sequent entrainment rate is no larger than a ~ 0.05. This 
reduced rate probably indicates that magnetic fields at the 
bubble/ICM interface suppress the growth of fiuid instabil- 
ities, while still permitting a small amount of thermal ma- 
terial to enter the bubble due to cross-field diffusion and 
magnetic reconnection events (e.g. |Pop3l2010l ). In any case, 
the large opening angles of the bubbles cannot be due to 
entrainment. Alternative possibilities include: 

(1) The bubbles are hydrodynamic in origin and the ap- 
parently large opening angle is partia l ly due to projection 
effects, as suggested bv lBriiggen et al.l l|2009l '): 

(2) The bubbles are not hydrodynamic in nature and 
are formed by a different infl ation proces s , e.g. current- 
dominated jets, as sugg ested bv lDiehl etall l|2008l ). 
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Figure 4. Plot showing tlie fraction of internal bubble pressure contributed by entrained material for different entrainmcnt rates: 
a = 0.15, 0.05, 0.01. By definition, the pressure fraction for the non-entraining bubble is always zero. 



6 DISCUSSION 

Having examined the three main mass transport processes 
associated with a buoyant bubble individually, we now con- 
sider the three processes jointly and ask: 

(1) how much material in total can actually be transported 
per bubble? 

(2) how frequently must bubbles be inflated in order for the 
mass outflow to balance the cooling inflow? 

Regarding the first question, we have already estab- 
lished that the bubbles cannot be carrying very much en- 
trained material. The drift can carry upward an amount of 
ambient material comparable to 50% of the mass initially 
displaced by the bubble, or approximately ~ 10* M©. How- 
ever, most of this material will not be displaced upwards by 
more than a bubble diameter and will likely fall back to- 
wards the cluster centre. In addition, for clusters with a flat 
entropy core, the overall fluid displacement might be negligi- 
ble. This then means that the drift, though likely responsible 
for observable features like the cool optical filaments, is not 
an important mechanism for affecting mass outfiow. 

This leaves t he wake. As noted previously, fiuidised bed 
experiments (e.g. lYandlioOsI ) show that wake typically con- 
tains a approximately a quarter of the mass initially dis- 
placed by the bubble. Assuming that laboratory value of q 
can be applied to the ICM, AGN-blown bubble wakes could 
easily transport gAfdis.o ^ 10* M©, see equation (j4|. 

Moreover, the mass in the wake can be carried out to 
large distances. As indicated by equation (|24|l , the maximum 
bubble height is a strong function of q. For q ~ 0.24 (and 
770 = 100), depending on the temperature profile, the bubble 
can rise to at least S^o, where zq is the height at which the 
bubble was infiated. For zo ~ 10 kpc, the bubble will rise to 
50 kpc. 

Let us consider an AGN that infiates a sequence of bub- 
bles. We can write the time-averaged mass outflow rate from 
the cluster centre due to a train of bubbles as 



where r is the average time between bubble inflation. Since 
this outflow counteracts the inflow due to cooling, cumu- 
latively bubble mass transport may allow for the relaxing 
of the stringent (but commonly adopted) requirement that 
AGN heating must balance radiative cooling. In fact, there 
is some debate as to whether AGN heating in fact balances 
cooling, especially in the most lumin ous clusters . Estim ates 
of the time-averaged heating rate bv lBest et al. I l|2007l ). for 
example, indicates that AGN heating does not balance cool- 
ing in the most luminous clusterslfl As we show below, mass 
outfiow can play an important role in limiting the mass 
build-up at the cluster centre in such instances. 

To illustrate this, let us consider the following: in the 
absence of heating, the mass infiow rate due to cooling is 
related to the radiative losses of the ICM according to 



M,, 



1) /J-mp 



TicM fceT 



ix. 



(51) 



where T is the average temperature of the ICM. Stud- 
ies of the X-ray spectra of the ICM in cluster cores in- 
dicate that the actual rate at which the g a s is cooling is 



Mcoo i ^0.2Mcooi.o (e g. iPeterson et all 120011 : iKaastra et al] 
I2OOII : iPeterson et"aLll2003l ). If the mass inflow rate is bal- 
anced by the mass transport rate by the bubbles: Mcooi ~ 

A/out. 



The actual gas cooling rate, Afcooi, 
difference between cooling and heating: 



is the result of a 



Mcool 



(TicM - 1) /imp 



(Lx - -Eheat) , (52) 



where Bhcat = ahcat-E'bub,o/-r. Here, £'bub,o is the enthalpy of 



Mout ~ qMnsfi/r, 



(50) 



^ We note that iDunn fc FabianI ll2008t) claim a close balance be- 
tween heating and cooling but in their comparison, they refer to 
the heating rate due to an individual bubble. Averaging over the 
time between successive bubbles will lower the power estimate, 
leading to a discrepancy be tween heating and cooling in agree- 
ment with lBest et al.l l[20o3). 
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an AGN-inflated bubble (see equation|3]) and cthoat allows for 
the possibility that the bubble enthalpy may un der-represent 
the total energy injected by the AGN (e.g. iNusser et al.l 
I2OO6I : iBinnev et al.lf2007h . especially if the bubbles are in- 
flated supersonically, which would engender shocks that en- 
hance the heating. 

Combining equations (jSOp — (|52p and taking Ficm = 
5/3, Tb = 4/3 and T/Tq ~ 3, we find that the build-up 
of excessive amounts of cold gas in the central galaxy can 
be prevented if 



and 



Mo, 



-Cihcat ~ JTi \ T-t^X 

[Iq + QhoatJ 

2q 



[2q + Qfhcat] 



A'/cool,0- 



(53) 



(54) 



If the AGN does no heating (i.e. aheat = 0), mass trans- 
port by the wake can, by itself, prevent rate of cold gas in 
the central cluster galaxy as long as the bubble recurrence 
timescale is 



1 X 10 



7 / -Bbub.o 
V 1059 erg ) V 10*4 ergs 



yrs, 



(55) 



where 10^® erg is the typical enthalpy of bubbles in a cluster 
with luminosity L x ~ 10*** ergs~^ (e.g. iBirzan et al]l2004l : 
iDunn et al.ll2005l ) and we have taken q « 0.24. Of course, 
the corresponding Mcooi will exceed the limits derived from 
the X-ray spectra. 

If, in addition to mass transport, the bubbles also heat 
the ICM such that the amount of energy deposited in the 
ICM is twice the enthalpy of the bubble (i.e. Qheat = 2), the 
rate of cold gas build-up can be prevented if the bubbles are 
inflated on a timescale 



5 X 10 



f -Ebub.O 

V 1059 erg ) V 10** ergs 



yrs. 



(56) 



It is worth noting that in this particular case, the time- 
averaged heating rate does not balance the cooling rate, 
Sheat ~ 0.8Lx, and the mass cooling rate is the maximum 
allowed, Mcooi ~ 0.2Mcooi,o, but mass build-up in the cen- 
tral cluster galaxy is successfully avoided by mass transport 
due to the bubbles. 

The values of r obtained from equation (|56|l are 
reasonably consistent with indirect observational esti- 
mates, though th ere are large uncertainties. For example, 
iBest et al.1 (|2007l ') found radio AGN duty cycles of ~ 30% 
in Brightest Cluster Galaxies (BCGs), which corresponds to 
an typical 'observed' time between outbursts of robs ~ 3ton, 
where ton is the duration of an AGN outburst. There- 
fore, since ton see ms to vary between 10^ — 10* yrs (e.g. 
iBirzan et alll2004h . it seems possible that AGN-blown bub- 
bles can significantly reduce the rate at which cold gas col- 
lects in the galaxy. 

Finally, we note that throughout this work we have only 
considered bubbles that remain intact. Bubbles which de- 
form into vortex rings behave somewhat differently. In ad- 
dition, the fluid circulation associated with the vortex ring 
exerts a forc e that can lift up addit ional material behind the 
bubble, (see iPavlovski et al.l 20071). Howe ver, with the pos- 
sible exception of M87 ( Owen et al.ll2000l ) structures resem- 
bling vortex rings have not been observed, so this additional 



mechanism of mass transport might be unimportant except 
in rare cases. 



7 SUMMARY 

The aim of this article is to provide a better understand- 
ing of the mass transport processes that are attributable 
to AGN-blown bubbles in clusters of galaxies. Recently 
observed fila ments of cool materi a l behind AGN-blown 
bubbles fe.g. [ Conselice et al.l I2OOII : ICrawford et al. I I-'IHI.-J : 
iHatch et al.l 2006[ ) provide a strong indication that bubble- 
induced mass transport is operating at least within 10s of 
kiloparsecs of the cluster centre. The filaments likely rep- 
resent only a fraction of the actual mass transported by 
the bubbles and some of this mass can be carried out to 
~ 50 — 100 kiloparsecs. 

Mass transport by AGN-inflated bubbles not only im- 
pacts bubble dynamics — including the rate at which the 
bubble expands as it rises — and the efficiency with which 
the bubble will heat its surroundings, it also removes the 
amount of cool gas available for star formation in the cen- 
tral galaxy, alleviating the stringent requirement that AGN 
heating balance radiative cooling. The system can tolerate 
a lower heating rat e, implying that the systems identified 
bv lBest et al.l (|2007l ). in which the average heating rates fall 
short of the cooling rates, may still be balanced. For all 
these reasons, a better understanding of mass transport by 
the bubbles is important. In this paper, we have focused on 
three main mass transport mechanisms: drift, wake trans- 
port and entrainment. 



• Drift: drift arises when an ascending bubble imparts 
an impulse on the ambient ICM ahead of it, causing a net 
upward displacement of the ambient material behind the 
bubble. In a gravitationally stratified environment, the ma- 
terial in the drift does not keep up with the bubble. It slows 
down and may even fall back to its original position. 

• Wake transport: a rising, deformable bubble also 
forms a cavity (wake) at the rear of the bubble. The wake is 
filled with some of the ambient medium displaced when the 
bubble first starts to rise, and this material will be trans- 
ported outward as the bubble rises. 

• Entrainment: during formation and possibly, during 
the rise phase, some of the ambient ICM may be incorpo- 
rated into the bubble body. This process is referred to as 
entrainment. The displacement of the bubble to larger radii 
in the cluster results in the transport of this entrained ma- 
terial. 

Our most important findings regarding the above three 
processes can be summarized as follows: 

(i) The upward displacement of the cool gas at the clus- 
ter centre in the form of drift offers a compelling explana- 
tion for the recently discovered cool filaments observed to 
be stretching between AGN-inflated bubbles and the cluster 
centre. The kinematic character of the filaments agrees well 
with the expected kinematic signature of drift. Our estimate 
of a drift mass of ~ 10* M© also agrees well with the esti- 
mate of mass in the cool filaments. We generically expect 
such filaments to form in clusters that either have entropy 
profiles that increase with radius (as is common in classical 
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cooling core clusters) or had such an entropy profile at the 
time when the bubble started its ascent. 

(ii) The wake is the most promising mechanism for trans- 
porting material out of the cluster centre and out to dis- 
tances of ~ 100 kpc, depending on where the bubble is 
inflated and the ICM's radial density profile. (In the ab- 
sence of the wake, the equilibrium radius of the very same 
bubble greatly exceeds the cluster's virial radius.) Our es- 
timates suggest that typically a bubble will lift ~ 10* M© 
in its wake, which significantly weighs down the bubble. A 
hot air balloon with a weighted basket is an apt metaphor. 
This additional mass not only limits the maximum radius 
to which a bubble can rise but also significantly reduces the 
energy that can be extracted from the bubble. In a classi- 
cal cool core cluster, a weighted bubble will deposit into the 
ICM only ~ 18% of the energy released associated with an 
unwei ghted bubble. 

(iii) iBriiggen et al] (|2009l ) recently suggested that en- 
trainment and projection effects may provide an explanation 
for the large apparent opening angles subtended by AGN- 
inflated bubbles. The presence of entrained material does 
increase the bubble's expansion rate. However, the entrain- 
ment rate required to account for the observations would 
also rapidly reduce the density contrast between the bubble 
and its surroundings. In a typical cluster, a bubble would 
only need to rise a distance equivalent to the initial bub- 
ble radius before it w ould become i ndisti nguishable from 
its surroundings. Since iBirzan et al.l (|2004l ) flnd numerous 
bubbles located more than 20 kpc from the cluster cen- 
tres, we conclude that the bubble does not entrain mate- 
rial efficiently. There are two related implications: firstly, 
the large opening angle must be due to some other pro- 
cess, such as cu r rent-d ominated AGN outflow model (e.g. 
iNakamura et al.l l2006|j_ as well as possib le projection ef- 
fects, as suggested bv lBriiggen et al.l (|2009l ). Secondly, some 
phenomenon prevents the formation of normal hydrody- 
namic instabilities and, therefore, suppresses entrainment. 
As a possible resolution, we note that bubble whose sur- 
face is draped by magnetic fields would behave as required. 
Whether such a bubble would be consistent with other con- 
straints is a subject for further explorations. 

(iv) With the introduction of mass transport of gas out 
of the cluster centre by the bubble, we flnd the clusters of 
galaxies can tolerate an average heati ng rate that is less than 
the cooling rate, which according to lBest et al.l l|2007h is a 
common occurrence in luminous clusters, and prevent gas 
build-up in the central galaxy. The outward mass transport 
rate is largely determined by the wake mass and by the av- 
erage time between bubble inflation. Typically, the bubble 
recurrence timescale is r ~ 10^ — 10* yrs, which is reason- 
ably consistent with observational estimates. 



through the Discovery Program and also expresses his grat- 
itude to J. Criswick for his generous support. AD thanks 
CITA for financial support through a National Fellowship 
and an NSERC grant to Don VandenBerg is also acknowl- 
edged. The authors wish to thank the anonymous referee for 
insightful comments that improved this work. 



8 ACKNOWLEDGEMENTS 

ECDP would like to thank the Department of Foreign Af- 
fairs and International Trade for funding through a Gov- 
ernment of Canada Post-Doctoral Research Fellowship and 
CITA for funding through a National Fellowship. He would 
also like to acknowledge informative discussions with Jim 
Hinton and Ian Fames. CP thanks the STFC for financial 
support. AB acknowledges research support from NSERC 



© 2010 RAS, MNRAS OOO.fTlfTSl 



16 E.C.D. Pope, A. Babul, G.Pavlovski, R.G. Bower, A. Dotter 



APPENDIX A: LIST OF PARAMETERS AND 
THEIR DEFINITIONS 



Parameter 


Definition 


-Ebub.O 


Till J-11 ■ J'J-1 -CJ- 

bub Die entiiaipy immediateiy atter 




inflation 


P 


/Till 1 T/^TV T 

pressure (bubble and iCM are m 




pressure equilibrium) 


Po 


pressure where bubble was inflated 


Vbub 


bubble volume 


Vbub.O 


initial bubble volume 


Pbub 


bubble density 


Pbub,0 


initial bubble density 


Tbub 


bubble temperature 




bubble adiabatic index 


Tb.o 


''I'll 111 1*1 I" • 1 

mitial bubble adiabatic index 


PiCM 


ICM density 


PiCM.O 


central ICM density 


po 


ICM density where the bubble was 




inflated 


TicM 


ICM temperature 


To 


ICM temperature where bubble was 




inflated 


Ti,T2, ZT 


parameters lor ht to ICM tempera- 




ture profile 


TlCM 


ICM adiabatic index 


Vo 


defined as po/pbub,o 




rp I" 1 ' 1 1 I r Till 

eftective density contrast tor bubble 




and wake system 


»7cff,o 


effective density contrast for bubble 




and wake system at the beginning of 




the rise phase 


a 
P 


exponent for /3-law ICM density 


Z 


radial coordinate of cluster 


2lCM,0 


scale height ot p-Iaw ICM density 


20 


height at which the bubble is in- 




flated 




1 ■ 1 J 1 1 ■ 1 ry J • 1 "J f* 

height at which effective density ol 




illlll / 1 J 111 

the bubble/wake systems equals the 




ICM density 


P(2eq) 


pressure at 2cq 


TlCM(2eq) 


ICM temperature at Zeq 


Vdrift 


drift volume 




wake volume 


Mbub 


bubble mass 


Mdis 


ICM mass displaced by bubble 


Mdis.o 


initial ICM mass displaced by bub- 




ble 


Mwake 


mass of bubble wake 




fraction of initial displaced mass in 




the wake 


k 


1' r ^ ' Ci 11 111 1 

ratio ol drift and bubbfe volumes 


9 


gravitational acceleration 


w 


bubble velocity 


^terminal 


terminal bubble velocity 


r 


bubble radius 


ro 


initial bubble radius 


flbub 


bubble surface area 


e 


turbulent energy dissipation rate 




per unit mass 


N 


Brunt- Vaisala frequency 
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Parameter Definition 
(cont...) 

a entropy index of ICM 

Cd drag coefficient for bubbie 

Lo Ozmidov iengtfi 

Lbuoy buoyancy lengtli 

L entrainment iengtli 

Cturb numerical turbulence constant 

Ctcrm drag coefficient-dependent 

Re Reynolds number 

n number of computational cells 

a entrainment coefficient 

Fi current bubble adiabatic index 

/i fraction of bubble pressure con- 

tributed by entrained material 
AV volume increment of entrained ma- 

terial 

AM mass increment of entrained mate- 

rial 

AVi volume entrained during ith 

timestep 

Vi bubble volume in ith timestep 

r effective adiabatic index 

Cp^bub specific heat at constant pressure of 

bubble 

CpjcM specific heat at constant pressure of 

ICM 

9 opening angle of bubble 

Lx X-ray luminosity of ICM 

T characteristic temperature of ICM 

Mcooi.o classical mass flow rate associated 

with Lx and T 
Mcooi the actual cooling flow 

Qbcat multiple of bubble enthalpy that 

represents total energy injected by 

AGN 

r average time between bubbles being 

inflated 

Tobs observationally-determined time be- 

tween bubbles being inflated 
ton duration of an AGN outburst 

fca Boltzman constant 

/imp mean mass per particle 

A cross-sectional area of bubble 
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